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ABSTRACT 

We report an X-ray spectral study of the transverse structure of the Centaurus A jet using new data from the 
Chandra Cen A Very Large Project. We find that the spectrum steepens with increasing distance from the jet 
axis, and that this steepening can be attributed to a change in the average spectrum of the knotty emission. Such 
a trend is unexpected if the knots are predominantly a surface feature residing in a shear layer between faster 
and slower flows. We suggest that the spectral steepening of the knot emission as a function of distance from 
the jet axis is due to knot migration, implying a component of transverse motion of knots within the flow. 

Subject headings: galaxies: active — galaxies: individual (Centaurus A, NGC 5128) — galaxies: jets — 
X-rays: galaxies 



1. INTRODUCTION 

In relatively few extragalactic radio jets can the kinematics 
of the flow on kpc scales be studied with any degree of cer- 
tainty. Most arguments are in direct, relying either on beam- 
ing statistics for samples (e.g.. [Hardcastle et al.ll2003al) or, for 
powerful jets oriented along the line of sight, on models re- 
lying on an as sumed inverse- Compton origin for the X-ray 
emission (e.g.. lSchwartz et aTl 2006). For the few low-power 
radio galaxies with heavily studied, straight, radio jets and 
counterjets, kinematic models hav e been constructed b ased on 
the jet-counterjet asymmetry (e.g.. lCanvin et al.ll2005h . These 
models have supported jet deceleration through mass entrain- 
ment, and ongoing work has already shown consistency be- 
tween the density and pressure model inferred from the jet 
and the properties of th e external X-ray-emitting gas for 3C 3 1 
(iLaing & Bridlell2002h . 

In models of mass entrainment, the outer parts ("sheath") 
are decelerated before the inner ("spine"), consistent with a 
range of observational evidence at radio frequencies (iLaingl 
[1996). AppHed to more central regions, the consequence that 
emission from a slower sheath becomes relatively more im- 
portant in jets at larger angle to the line of sight then re- 
solves difficulties in models t hat unify BL Lac obje cts with 
low power radio galaxies (e.g.. lChiaberge et"ani2000h . 

Chandra has not only allowed confirmation t hat X-ray syn- 
chrotr on emission in low-power jets is common (IWorrall et al.l 
1200 Ih . but has also resolved transverse structure in the near- 
est and brigh test sources. A part icularly interesting example 
is NGC 315 (IWorrall et al.ll2007b . Here the diffuse emission 
contains a knotty structure in the radio and X-ray that ap- 
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pears to describe an oscillatory filament. Although the struc- 
ture could be the result of a chance superposition of non- 
axisymmetric knots, the level of coherence led us to suggest 
that the knots might be predominantly a surface feature resid- 
ing in the shear layer between the spine and sheath. If this lat- 
ter interpretation is correct, we might expect the X-ray spectra 
of the knots to be similar across the transverse width of the jet. 

In the case of NGC 315 the distinct knotty emission is only 
about 10% of the total emission in X-rays and radio along the 
^ 2.5 kpc of projected jet length over which it is detected, 
and with a source distance of ~ 70 Mpc the observations did 
not allow us to measure separate X-ray spectra for the knots 
and diffuse emission. At 3.7 Mpc** (1 arcmin is 1.1 kpc 
projected length), Centaurus A is a much closer example of 
a low-power radio galaxy whose inner jet shows bright re- 
solved X-ray knots and diffuse emission over a similar pro- 
jected length scale to NGC 315 (e.g. Ki-aft et al. 2002). Here 
we report an X-ray spectral study of the transverse structure 
of the jet, and in particular we investigate whether the Cen A 
knots are likely to occupy a shear layer between faster inner 
and slower outer flows. 

2. OBSERVATIONS & ANALYSIS 

The Chandra data used here are from the Cen A Very Large 
Project (Cen A-VLP), comprising six deep (~ 100 ks) expo- 
sures (OBSIDs 7797, 7798, 7799, 7800, 8489, 8490) of the 
front-illuminated CCD chips of the Advanced CCD Imag- 
ing Spectrometer (ACIS) made between 2007 March 7 and 
May 30. Standard methods'* have been followed for data re- 
processing and grade selection, and small astrometric correc- 
tions were made to register the six e xposures to a common 
frame using detected point sources (iHardcastle et all 120071: 
Jordan et al. 2007). The 0.8-3 keV merged image of the in- 
ner ~ 4' of the X-ray jet is shown in Figure[T] 

We are interested here in the jet after it has become radia- 
tive ly bright in both radio and X-ray (more than 21" from 
the core , and after knots AXl A, AXlC in the terminology of 
IHardca stle et al. 2003b) until its major decline relative to the 
radio (at 141"). Motivated by NGC 315, where the fast 
spine is approximately 60% of the jet width, we have defined 

* The average of five distance indicators, see §6 in lFerrarese et all 120071) 
' http://asc.harvard.edu/ciao/ 
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Fig. 2. — Image from Fig. [T] with regions superposed, and rotated to 
the sky orientation. The regions, called the spine, sheath and background, 
working outwards from the jet axis, are pie slices from a base position of 
R.A. = 13'^25™26?98, Dec. = -43° 01' 14'.' 06 (not the core) running from 
position angle (pa) 51?2-57°8 for the spine, pa 49°-51?2 and 57? 8-60° for 
the sheath, and pa 31°^7° and 63°-79° for the background with the exclu- 
sion regions shown. From the core the inner spine and sheath regions are at 
a distance of 21" to 66", and the outer region runs from 66" to 141". We 
call the region marked in blue, near the outer extent of the inner region, the 
"bright" region. 

regions shown in Figure |2] While not wishing to suggest a 
priori that the Cen A jet does have a slower sheath and faster 
spine, for shorthand we call these regions the spine, sheath 
and background, working outwards from the jet axis. Com- 
pact sources excluded from the background regions are also 
marked. A division of the spine and sheath into two regions 
with length down the jet is shown. We refer to these as the 
inner and outer regions, and the boundary was chosen so that 
there is a comparable number of total counts in each region. 

The jet is positioned on a single ACIS I chip in all six 
exposures. The boundary between the outer and inner re- 
gions lies between 1.4 and 4.0 arcmin off axis depending 
on OBSID, giving a point-spread-function (PSF) half-power- 
diameter (HPD) at this location of better than about 1 arcsec 
at the ^ 1.5-keV energy of maximum count rate per keV. The 
readout streak from the Cen A core runs close to the SE sheath 
region in OBSID 7797, and so we exclude this observation in 
spectral fitting to the SE sheath and inner spine. We extract 
data and responses for identical sky regions in the different 
OBSIDs, and we fit the same model simultaneously to the 
data sets. Spectra are binned to a minimum of 30 counts per 
bin so that the statistic can be applied. 

In this work we are searching for statistical trends in spec- 
tral index in knotty and diffuse emission transverse to the jet 



axis. Knots are identified as bright enhancements, but we be- 
lieve them to be composed of multiple distinct physical struc- 
tures, each potentially with a different emission spectrum. 
The precise definition of knots is therefore not cruci al, and 
we hav e ad opted the previ o usly p ublished lists of Kraf t et all 
(120021) and iKataoka et al.1 (l2006l) . which combined give 31 
knots within our total spine and sheath regions. The knots 
were examined visually, and small centroid shifts (within un- 
certainties) were made to give best alignment in the new data. 
A circle of nominal radius 2.4 arcsec was used to mask each 
knot. Four additional weak but compact enhancements seen 
in the new data were also masked as knots with circles of 
radii 1.2 arcsec. Because of our interest in statistical trends, 
for knots that overlap regions we divide the counts accord- 
ing to where they were detected. When extracting spectra of 
the knots the diffuse (residual) emission in the correspond- 
ing region is used as background. The outer sectors shown in 
Figure |2]provide background for the diffuse emission. 

All our spectral models are single-component power laws 
with free hydrogen column density, A^h, over the energy 
range 0.4-10 keV. Intrinsic A^h in Cen A over and above the 
Galactic value of 8.4 x lO^o cm'^ dPickev & Lockmanll990l) 
causes reductions in the count rate in soft X-ray images, con- 
sistent with the dust lanes that cross the jet roughly transverse 
to its axis. Since our regions extend hundreds of pc along 
the jet, we are combining data with different intrinsic A^h- 
This is expected to give a bias towards finding slightly flat- 
ter power-law slopes than the true slope if the jet emission is 
characterized by a single spectrum over its length. However, 
the gradients in A'h are principally down the jet and not across 
the jet, and so provide similar weightings when we compare 
sections transverse to the jet axis. We can thus use the fit- 
ted A'h and spectral index for our transverse sections of the 
jet as proxies of systematic trends in the spectra. Appreciable 
cross-jet differences could then allow us to conclude that there 
are real variations in the X-ray spectrum across the jet. Self- 
consistency tests (see verify this assertion. We also find, 
through comparison with the spectra of indi vidual knots and 
sections of diffuse emission along the jet ( Hardcastle et al.l 
l2007h . that any bias in spectral slope in our analysis is small. 

In the following, spectral index a (one less than X-ray pho- 
ton index) is defined in the sense that flux density is propor- 
tional to i^^". 

3. RESULTS 

Our spectral results are summarized in Table [T] and we re- 
fer to them here by row number The quoted uncertainties are 
90% for two interesting parameters, so where pairs of rows 
disagree in both parameters there is at most a 1 % probability 
of agreement by chance. Uncertainties in a and A'h are cor- 
related, and so the contours must be examined to assess 
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TABLE 1 

Chandra FITS TO A SINGLE-COMPONENT POWER LAW 



Row 


Component/ 


region % counts 


region % area 




A^H (10^1 cm ^) 


X^in/dof 


1 


Spine 


100 


100 


n 00+0.05 


1 54+0.18 


961/895 


2 


Sheath 


100 


100 


0.97in? 


1.40 ±0.28 


527/512 


3 


NW sheath 


100 


100 




1 47+0.46 
r.^/ —0 44 


318/277 


4 


SE sheath 


100 


100 


0.99 ±0.11 


1 51+0.45 


295/282 


5 


Spine knots 


64 


25 


0.60 ±0.06 


77+0.29 
—0 28 


602/657 


6 


Sheath knots 


49 


17 


0.73l«;i? 


1 19+0^53 
— 1 06 


244/254 


7 


Spine diffuse 


36 


75 


1.24 ±0.10 




471/407 


8 


Sheath diffuse 


51 


83 


1.20 ±0.14 


1 09+0.51 

l-<5''_0 50 


338/308 


9 


Inner spine knots 


77 


47 


0.63 ±0.08 


.^"^^—O 42 


474/520 


10 


Inner sheath knots 


61 


33 




1 64+°-*0 


221/185 


11 


Outer spine knots 


34 


18 






209/207 


12 


Outer sheath knots 


27 


13 


o.8oig:i° 




60/71 


13 


Inner spine diffuse 


23 


53 


0.96t°:}? 


9 05+0.83 


182/168 


14 


Inner sheath diffuse 


39 


67 


u.au_Q 22 


9 — n + 1.27 

^- '■^-1.16 


136/121 


15 


Outer spine diffuse 


66 


82 


1 qc+0.12 




318/336 


16 


Outer sheath diiTuse 


73 


87 


1 qo+0.19 


1-^1-0.59 


191/191 


17 


"Bright" spine knots 


88 


57 


•^•'^-o.io 




229/249 


18 


"Bright" sheath knots 


82 


60 




1 19+0.96 
l.l^-n.90 


147/124 



' En'ors are 90% for 2 interesting parameters 




1.6 



1.8 



2.2 




Photon Index (a + 1) 



Fig. 3. — Confidence contours (Icr, 90% and 99% for two interesting pa- 
rameters) showing a significant difference between the spine inner knots (left 
contours) and sheath inner knots (right contours): rows 9 & 10 of Table [T] 
The y axis begins at Galactic A'^jj- 



fully the spectral disagreement. We show such plots only for 
the comparisons of most interest (Figs.[3]&|4l). 

We see (Table [T] rows 1 & 2) that the sheath spectrum is 
steeper than the spine. Further investigation allows us to find 
the origin of this result. The NW and SE regions of the sheath 
(rows 3 & 4) give excellent agreement in a and iVn, indicating 
that it is unlikely that the result is from isolated features, and 
confirming the absence of strong gradients in A^h transverse to 
the jet. Separating the knot and diffuse emission in the spine 
and sheath (rows 5 to 8) shows that while the knots have flatter 
spectra than the diffuse emission, the steeper spectrum in the 
sheath is caused solely by differences in the knots. In the two- 
dimensional plane of a and iVn, rows 5 and 6 disagree. 

Having established that the knots in the sheath have a 
steeper spectrum than in the spine, we checked the inner and 
outer regions to see if there is any trend with distance from 
the core. Comparing rows 9 and 10 with 11 and 12, we see 
the expected preference for a higher A'h in the inner regions, 
but more interestingly we see that it is only for the knots 
in the inner regions that the spine and sheath spectra differ 
(Fig-E). The diffuse emission, in contrast, shows remarkable 
spectral agreement between the sheath and spine in the inner 
region (rows 13 & 14). To check that the contour separation 



Photon Index (a + 1) 

Fig. 4. — Same as Fig.|3]for the bright region spine knots (left contours) 
and sheath knots (right contours): rows 17 & 18 of Tabled] 

in Figure|3]is causal rather than statistical we have drawn ran- 
dom sets of the required size from the modelled underlying 
spectral-index distribution deduced from all knots in the inner 
region. The average spectral index has an error on the mean of 
±0.07 for knots matching in number those in the inner spine, 
and ±0.11 for the inner sheath. The observed spectral-index 
difference, Aa — 0.31 (Fig. |3]l, therefore significantly ex- 
ceeds the difference expected from random subsamples of the 
distribution. 

Since the spine inner jet knots tend to lie in regions of 
higher A'^h than the sheath inner knots (Figs|2]and[3l), we have 
repeated the analysis only for the "bright" region, marked in 
blue on Figure |2] This region i s beyond the main dust lanes 
that cross the jet (e.g., fig. 1 of iHardcastle et alll2007l) . Re- 
sults are consistent with sheath knots having steeper spectra 
than spine knots (Fig. |4] and rows 17 & 18). Simulations 
show that the bias (§2) on the absolute value of mean spectral 
index from the A'h distribution over the inner jet is ~ 0.05. 

We investigated the possibility that the spectral flattening 
in the inner spine might result from the inner knots being 
brighter, if brightness is associated with particle-acceleration 
efficiency. Figure |5] shows the knot surface brightness and 
hardness across the inner jet. The bin-width is 1 degree in 
position angle, which corresponds roughly to the HPD at an 
average jet-axis distance, and so the points are not highly cor- 
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Fig . 5 . — Surface-brightness profiles with position angle of the knot emis- 
sion from the inner jet in a hard energy band (upper panels) with the hardness 
ratio below. The X-ray emission is softer towards the jet edges (x^ = 116 
for 1 1 dof for a fit to constant hardness). There is no trend for regions of 
higher integrated knot suii'ace brightness to show harder spectra. 

related. The hardness ratio shows the decrease with increasing 
distance from the jet axis that is anticipated from our spectral 
results, with a high degree of confidence that the hardness is 
not constant. The gradual decrease of hardness with distance 
from the jet axis shows that our finding of a spectral differ- 
ence in the spine and sheath regions was not dictated by the a 
priori definition of the region boundaries. 

Perhaps surprisingly, a comparison of the panels of Figure|5] 
shows that there is no trend for regions of higher integrated 
knot surface brightness to be harder in spectrum. We have 
made plots similar to Figure|5]for the inner and outer diffuse 
emission, and the outer knot emission. The surface-brightness 
profiles of diffuse emission show some lumpiness, suggest- 
ing that some of what we call diffuse emission is composed 
of lower-surface-brightness knots, but all three hardness-ratio 
profiles are consistent with being constant across the jet (x^ 
values of 10.5, 7.9, and 14.8 for 1 1 degrees of freedom (dof), 
as compared with =116 in Fig.|5]l. 

4. CONCLUSIONS 

If the knots in Cen A were predominantly a surface feature 
in a shear layer between the spine and sheath, as we suggested 
might be the case for NGC 315, one would expect the knot X- 
ray spectra to be similar across the jet. This is not what is 
seen. In the inner jet, between 21" and 66" from the core, we 
have observed in a statistical sense that the X-ray knot spectra 
steepen farther from the jet axis. No such trend is seen in the 
more diffuse X-ray emission. The results suggest that knots 
are distributed throughout the jet volume, with a trend towards 
steeper spectra in the outer regions. The trend will then be 
underestimated in our results because of projection effects. 



Although the spectral trends disfavor knot confinement in 
a shear layer, they imply that the electron energy distribu- 
tions in knots towards the center and edge of the jet are dif- 
ferent. Electrons with sufficient energy to emit synchrotron 
X-rays must be accelerated locally to the knots, since their 
energy loss timescale is short compared with any dynamical 
timescale associated with fluid motion. Thus the conditions 
for particle acceleration differ across the jet. The knots might 
arise from intruders from outside travelling at speeds close to 
the stellar velocity dispersion, or internal flow irregularities 
which are most likely to arise in the central part of the flow 
and might be expected to move quickly towards the edge of 
the jet under the influence of circulation in the shear layer. 
In either case a transverse velocity structure might cause the 
strength of shocks near the knots to vary across the jet, and 
since the kinetic energy density of the flow should be higher 
nearer the jet axis, we might also expect velocity irregularities 
to generate stronger turbulent cascades near the axis than near 
the jet edge, with con sequent efficient particle acceleration 
(iManolakou et al.lfl999l) and a flatter on-axis X-ray spectrum. 

As the jet expands, the lower on-axis flow velocity, as sug- 
gested by the less collimated appearance of the jet further 
from the core, would then imply less X-ray spectral differ- 
ence between sheath and spine, consistent with our results. 
Indeed, a spectral steepening between the inner and outer re- 
gions in diffuse X-ray emission may relate to such a decline 
in flow speed along the jet (rows 13 to 16 of Table[T]). 

It is interesting to ask whether it will be possible to de- 
tect knot migration in Cen A. Expansion of the jet along 
its axis, with apparent motions of about 0.5c, has been mea- 
sured in radio datafrom the Very Large Array over 10 years 
dHardcastle et al.l l2()03bt) . but lateral motion is more diffi- 
cult to measure. It has been detected in the radio in M87 
(iBiretta et al.lll995h . We speculate that transverse motions 
should be present in the knots and are potentially detectable 
if knots are flow irregularities, although future observations 
may be required. 

These new observations of Cen A have permitted the first 
X-ray spectral study of a radio jet in the transverse direction, 
separating knotty and more diffuse emission. We have found 
the striking result that where the X-ray emission is brightest, 
at ^ 0.4 to 1.2 kpc from the core, the spectrum in the knotty 
emission is steeper at larger distances from the jet axis. We 
speculate that this is the result of lateral knot migration. 
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